Complex Spontaneous Flows and Concentration Banding in Active Polar Films 
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We study the dynamical properties of active polar liquid crystalline films. Like active nematic 
films, active polar films undergo a dynamical transitions to spontaneously flowing steady- states. 
Spontaneous flow in polar fluids is, however, always accompanied by strong concentration inhomo- 
geneities or "banding" not seen in nematics. In addition, a spectacular property unique to polar 
active films is their ability to generate spontaneously oscillating and banded flows even at low 
activity. The oscillatory flows become increasingly complicated for strong polarity. 



Active materials are a new class of soft materials main- 
tained out of equilibrium by internal energy sources. 
There are many examples in biological contexts, includ- 
ing bacterial colonies [l| , purified extracts of cytoskeletal 
filaments and motor proteins [2], and the cell cytoskele- 
ton A non-biological example is a layer of vibrated 
granular rods 0] . The key property that distinguishes ac- 
tive matter from more familiar non-equilibrium systems, 
such as a fluid under shear, is that the energy input that 
maintains the system out of equilibrium comes from each 
constituent, rather than the boundaries. Each active par- 
ticle consumes and dissipates energy going through a cy- 
cle that fuels internal changes, generally leading to mo- 
tion. The experimental systems studied to date typically 
consists of elongated active particles of two types: polar 
particles, with a head and a tail, and apolar ones that 
are head-tail symmetric. Active suspensions can then 
exist in various liquid crystalline states, with novel struc- 
tural and rheological properties 0, 0] . Apolar particles 
can form phases with nematic order, characterized by a 
macroscopic axis of mean orientation identify by a unit 
vector n and global symmetry for n — > — n, as in equilib- 
rium nematic liquid crystals. Polar particles can order in 
both nematic and polar phases. The polar phase is again 
characterized by a mean orientation axis p, but p / -p. 
The protein filaments which are the major component of 
cell extracts are generally polar and these extracts can 
therefore have both nematic and polar phases. 

Conventional liquid crystals exhibit a rich non- 
equilibrium behavior when subject to external forcing, 
such as shear or applied magnetic and electric fields. This 
includes transitions to stable statically distorted defor- 
mations of the director field (Freedericksz transition Q), 
shear banding [8], and even the onset of turbulent and 
chaotic behavior in the presence of shear [9j]. Active liq- 
uid crystals exhibit a similar wealth of phenomena due 
to internal forcing, i.e., spontaneously. A striking prop- 
erty of active nematic liquid crystal films is the onset 
of spontaneous flow above a critical film thickness first 
identified by Voituriez et al [lOj]. This phenomenon is 
analogous to the Freedericksz transition of a passive ne- 
matic in an applied magnetic field, but the flowing state 



is driven by the internal activity of the system - hence the 
name. This prediction was obtained by analytical stud- 
ies of the phenomenological hydrodynamic equations of 
an an active nematic film in a one-dimensional geome- 
try. More recently, Marenduzzo et al [ll| have employed 
hybrid lattice-Boltzmann simulations to study the active 
nematic hydrodynamics in both ID and 2D geometry and 
have mapped out the phenomenon in parameter space. 

In this letter we show that active polar fluids exhibit 
an even richer behavior. First, like active nematics, po- 
larized active liquid crystals exhibit steady spontaneous 
flow. Unlike active nematics, however, where the filament 
concentration remains practically uniform in the sponta- 
neously flowing state, spontaneous flow in polar fluids 
is accompanied by "concentration banding", i.e., a sharp 
gradient in the concentration of filaments across the film. 
The concentration banding is due to active couplings of 
concentration gradients to the polar director in the hy- 
drodynamic equations that are allowed only in fluids with 
polar symmetry. Upon increasing the magnitude of these 
polar couplings, the steady state becomes unstable and 
the system undergoes a further transition to a dynamic 
state with bands of oscillating concentration and orienta- 
tion. In the oscillatory regime, travelling bands nucleate 
and oscillate from one end of the film to the other. For 
even larger couplings the oscillatory behavior becomes 
increasingly complex with the appearance of multiple fre- 
quencies with incommensurate ratios between the periods 
of the orient ational and concentration oscillations. 

Hydrodynamic equations for a two component ac- 
tive suspension have been written down phenomenologi- 
cally [12] and derived from a microscopic model [l4]. The 
relevant hydrodynamic variables are the concentration c 
of filaments and the total density p and momentum den- 
sity g = pv of the suspension, with v the flow velocity. 
We consider an incompressible film, with p = constant 
and V • v = 0, and macroscopic dimensionless polarity 
P = |P|p, with direction characterized by a unit vector 
p = (cos sin the polar director. 

The hydrodynamic equations for two-component po- 
lar suspensions have been derived elsewhere by coarse- 
graining a microscopic model of the dynamics of inter- 
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acting motors and filaments . The active suspension 
is an intrinsically non-equilibrium system and cannot be 
described by a free energy. However, for clarity of pre- 
sentation we introduce here the equations phenomeno- 
logically and write all equilibrium-like terms (i.e., those 
terms that are also present in an equilibrium polar sus- 
pension) in terms of derivatives of a non-equilibrium ana- 
logue of a free energy, given by 
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with C the compressional modulus and K\ and K% the 
splay and bend elastic constants, both taken equal to K 
below. The last three terms on the right hand side of 
Eq. (pQ) couple concentration and splay and are present 
in equilibrium polar suspensions [B\ = B2 = B3 = B 
below). The dynamics of the concentration and of the 
polar director is described by 



d t c 



c(v - (3'c£ 2 V) + T'h 



d t + (v + Pc£ 2 P) • V 



LdijPj 



Xuij 



Pj + rhi-r'di^), (2b) 



, (2a) 



where I is the length of the filaments, uij = \{diVj + 
djVi) and Uij = \{diVj — djVi) are the rate-of-strain and 
vorticity tensors, and h = — 8F/SP the molecular field. 
Here T, V and V" are kinetic coefficients and (3 and f3' are 
active parameters. The equation describing momentum 
conservation is written in the Stokes approximation as 
0. The stress tensor cr 71 - is the sum of reversible, 
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dissipative and active contributions, 
The reversible part is written in an equilibrium-like form, 



with n the pressure and A and A alignment parameters. 
The dissipative part of the stress tensor is written as 
afj = rjUij assuming a single viscosity 77. Finally, there 
are additional stresses induced by activity given by 
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with n a the active part of the pressure. There are two 
contributions to the active stress tensor. The first (~ a) 
describes active stresses that arise from contractile (if 
a > 0) forces induced by activity. This term is present 
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FIG. 1: (Color online) (a) Sketch of the film geometry, (b-d) 
Solutions of Eqs. d3a}, J3bj) for A = 0.1, f = 0.3,^0 = 1, 
C = 0.5, T]F = 0.5, w = 0.13, a = 0.08 and variable p. 



in both nematic and polar liquid crystals and its effects 
have been studied before. The second term (~ f3") arises 
from "self-propulsion" of the active units and is exclusive 
to polar systems. The same mechanism is also responsi- 
ble for the "connective" terms proportional to (3 and f3' in 
Eqs. ([2a]) and (|2bp . For a motor /filament mixture, all ac- 
tive contributions are proportional to the mean rate Afi 
of ATP consumption, which is the internal driving force 
for the system, ^-type terms have dimensions of veloc- 
ity and have been estimated in the microscopic model as 
(3 ~ mwo, with fh a dimensionless concentration of motor 
clusters and uo ~ A/i the mean velocity at which motor 
clusters step along the filaments [l3[. Previous work on 
confined active films [Til ] has been limited to active 
nematics with all the f3 terms equal to zero. Here for 
the first time we incorporate the polar active terms and 
analyze their role in controlling non-equilibrium effects 
in active films. Hereafter we will assume (3' = (3" = (3. 

We consider a two-dimensional active polar suspension 
with polarization of uniform magnitude and discuss the 
dynamics of the hydrodynamic fields c, v and p. For sim- 
plicity we set |P| = 1. The film sits on a solid plane at 
y = and is bound by a free surface at y = L (Fig. [T^l) . 
The discussion below is easily extended to other bound- 
ary conditions. The film extends to infinity in the x di- 
rection and we assume translational invar iance along x. 
The Stokes equation requires d y cF yy = and d y a xy = 0. 
The first of these two conditions fixes the pressure in the 
film. The second, together with the boundary condition 
cr xy (L) = 0, requires a xy = constant = throughout the 
film. We also assume no-slip boundary conditions at the 
substrate, so that v x (0) = 0. 



It is convenient to work with dimensionless quantities 
by introducing the time scale ro = £ 2 /(TK). Letting z = 
y/£, r = t/ro, (j) = c/cq and specializing Eqs. (f2a|) and 
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FIG. 2: (Color online) On the left discrete Fourier transforms 
of 0{zl = L/2,t) (center of the film) for A = 0.1^ f = 0.3, 
D = 1, C = 0.3, nT = 0.5, w = 0.13, a = 0.1 and (3 = 15 (a) 
18 (b) and 20 (c). On the right space-time plots of <j>(z,r). 



(|2b|) to our quasi-one-dimensional geometry, we obtain 

<9 r = d z ^f3(j) 2 sin 6 + AS sin sin 2$ 

+ [L>(1 - £ sin 2 (9) - w cos 2 (9] d^j , (3a) 

d r = (l-w cos 2 0)<9 2 + \w sin 20(d z 0) 2 

- # sin 0<9 2 + w cos 0d z (j) - 2(1 - A cos 2(9) , (3b) 

where /? = /?c fro, w = 21B/K, £ = r /2 /(rr // ), and 
D = r^C/^oFif). The dimensionless rate-of-strain u = 
1u xy To can be obtained from the condition a xy = 0, 
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where 5 = ac 2 £ 4 /(TK) and C = fT'C/fTK) 2 . The hy- 
drodynamic equations for a two component active ne- 
matic suspension are obtained from the above by letting 
(3 = and w = 0. The terms proportional to w are 
also present in passive polar fluids as they arise from the 
fact that the polar symmetry allows the coupling propor- 
tional to B between splay and density fluctuations. The 
terms proportional to (3 are intrinsically non-equilibrium 
polar terms. Finally, the case of an incompressible one- 
component nematic fluid, investigated by Voituriez et 
al [10[ and by Marenduzzo et al [11], can be recovered 
from our equations by setting (3 = w = and assuming 
a constant concentration 6. 



Eqs. (|3bj) and (|3a|) are integrated numerically with 
boundary conditions 0(0, r) = 6(L/£,r) = 0, d z (j)(0, r) = 
d z <t>{L/i,r) = (i.e. j y (0,t) = j y (L,t) = 0). The ini- 
tial conditions on and (j) are chosen as random, with 
the constraint (0(^,0)) = and (</>(z,0)) = 1 where (•) 
stands for a spatial average. 

Steady spontaneous flow: Both the polar and the apo- 
lar systems exhibit a Freedericksz-like transition between 
a state where the director field is constant and parallel 
to the walls throughout the channel to a non-uniformly 
oriented state in which the system spontaneously flows in 
the ^-direction. The transition can be tuned by chang- 
ing either the film thickness or the activity parameter a. 
Fig [T] shows a numerical solution of Eqs. ([3aj) and (|3b|) 
for fixed 5 and variable /3, with L/i = 10. As the active 
velocity (3 is increased, the maximum tilt m decreases 
and the alignment is progressively restored. Remarkably 
the variation in the concentration <\> across the film is 
significantly stronger than in the apolar case (solid green 
curve in Fig. [T]) with a relative difference between the 
highest and the lowest values up to 50%. This "concen- 
tration banding" is a characteristic of polar active sys- 
tems. It is a consequence of the active (3' coupling in Eq. 
([2a|) resulting from self-propelled 'convection' of the ac- 
tive elements along the local polarization direction. The 
varying local polarization angle required for spontaneous 
flow therefore leads to an even stronger variation jn the 
local concentration. Close to the transition, 5 c (/3), the 
coupling between the polar director and concentration 
also leads to an asymmetric director profile across the 
film. We also point out that in the absence of this polar 
active term, there are equilibrium-like gradient couplings 
between local director and density. These are, however, 
much weaker since they occur at higher order in gradi- 
ents. In contrast, the active nematic shows a negligible 
concentration gradient even for anomalously large values 
of the contractile activity parameter, 5 ^> a c ((3). 

Spontaneous oscillations: Upon further increasing /?, 
spontaneous oscillations of <j){z : r) and 6(z,r) are ob- 
tained. The coupled dynamics of the two fields gives 
rise to travelling waves of concentration and orientation 
bands. Initially only one frequency is observed, but the 
oscillations become more complicated as (3 increases. Us- 
ing Fourier decomposition we find that this is due to the 
appearance of additional frequencies at different values 
of (3 for concentration and orientation bands (see Fig [2]). 
A phase diagram in the (a,/?)— plane is displayed in 
Fig. 03 It shows transitions between stationary (S) flow, 
spontaneous steady flow (SF) and spontaneous periodic 
(oscillatory) flow (PF). 

We can understand the phase behavior close to a 
stationary homogeneous state (0 = 0o, = 0, u = 
0) by expanding 0(z) = 56+(z)e iujt + 56-(z)e- iujt and 
(j){z) = 0o + S(f)+(z)e iu;t + 8(j)-(z)e- iujt . The boundary 
conditions require 4>±{z) = Y^=i a n cos (niriz/L) and 
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FIG. 3: (Color online) Phase-diagram in the (a, (3)— plane. 
The points have been obtained numerically using A = 0.1, 
f 0.3, D = 1, C = 0.3, f]F 0.5, iu 0.13 and L/£ 
10. For /3 < Ptp upon increasing a the system undergoes 
a transition between stationary homogeneous state (S) and 
inhomogeneous steady flow (SF). Above the "tricritical point" 
(TP) the spontaneous flow becomes oscillatory (PF). 



0±( z ) = S^Li sm (nir£z/L). The first unstable modes 
are af : bf and a linear stability analysis shows that there 
is a steady-state (u = 0) instability at 
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to a steady spontaneously flowing state with concentra- 
tion banding. Oscillatory modes with frequency uo c ~ 
</> (£/L)(uj0) 1 / 2 appear beyond a "tricritical point" (TP) 
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This linearized analysis predicts the positions of the 
"phase boundaries" S-SF and S-PF in quantitative agree- 
ment with the numerical solution. The boundary SF-PF 
has been obtained only numerically. The appearance of 
spontaneous oscillations results from the coupled motion 
of concentration and director orientation bands due to 
both the convective active polar coupling (/?) and the 
passive polar coupling (w) of director and concentration. 
Upon increasing the oscillatory behavior becomes in- 
creasingly complex, but we have not been able to observe 
fully fledged chaos for reasonable values of 0. 

We have studied the dynamical properties of thin films 
of active polar fluids and found a rich variety of com- 
plex behaviors which should be observable experimen- 
tally in polar active systems. Using microscopic models 
of motor-filament coupling, it was estimated in [13] that 
(3 ^> a in microtubules/kynesin mixtures, while (3 < a in 
actomyosin systems. This suggests that in- vitro micro- 
tubules/kynesin mixtures may be the best candidate for 



the observation of the oscillating bands predicted here. 
It should, however, be noted that filament treadmilling 
also leads to terms with polar symmetry at the contin- 
uum level, where in this case would be proportional 
to the polymerization rate [la]. The intriguing possi- 
bility that our findings may be relevant to treadmilling 
acto-myosin systems and therefore have implications for 
lamellipodium dynamics will be explored elsewhere. 
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